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Abstract 



A test mass, M, moving through an ambient medium of light particles with lower 
average kinetic energy than itself suffers a deceleration caused by its scattering of 
the light particles. The phenomenon is usually referred to as dynamical friction. 
The velocity, v, of the test mass decays on a timescale independent of v in the 
non-relativistic case. We derive expressions for dynamical friction in the case that 
the test mass and the light particles are relativistic, and that the scattering is 
weak (with impact parameter, b ^> M). In the case that the light particles are 
ultra-relativistic, and isotropic in the frame in which M moves with velocity v, we 
find an explicit expression for the dynamical friction. The well known factor of 2 
correcting the Newtonian scattering of photons to give the Einstein angle, 4M/6, 
has the largest effect on the resulting friction, which is modified by a factor of 
roughly 16/3j v over the simple non-relativistic case. In the non-relativistic case, 
the largest contribution to the friction comes from light particles moving slower 
than v. We find that this is not the case for ultra-relativistic scattering, essentially 
because the scattering angle is independent of v. Some astrophysical implications 
are discussed. 



1 Introduction 



Consider a test mass, M, moving through a medium of light particles. It will be 
useful to define the 'lab frame' as the one in which M moves with velocity v. In 
what follows we shall adopt natural units in which G = c = 1. Chandrasekhar 
(1943) derived a formula for the dynamical friction on M, which may be written 

dv 9 , C f(w)w 2 dw 

— = -lQ7r 2 AMm Jo V ; v, (1.1) 

dt v s 

where A is the Coulomb logarithm, and f{w) is the distribution function of light 
particles, assumed to be isotropic. In the case that v is much less than the average 
value of w, and / is Maxwellian with dispersion w, 
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— ~ -4vryl^-v. (1.2) 
at w A 

One important feature of (1.1) is that only the light particles with w < v contribute 
to the friction, to order A. If all the light particles were massless, there would be 
none with w < v } so a naive application of (1.1) would predict no friction at all 
in this order. Of course, we would not be surprised to find contributions at other 
orders, but in the relativistic case A can be very large, so it will be important 
to investigate the precise form of the dynamical friction. We shall show in the 
next section that the friction is actually enhanced if the light particles are ultra- 
relativistic. Our derivation is closer to that of Binney and Tremaine (1987) than 
to the original work by Chandrasekhar. 



2 Dynamical Friction 



Consider an event in which a light particle is scattered through a small angle, x, by 
gravitational interaction with M. The velocity change imparted to M in a single 
collision is of order m/M times that imparted to m. In the absence of a solution to 
the full relativistic two-body problem, we only know how to solve for the motion 
of m in the limit that this quantity is small, m/M <C 1. The motion of M will 
be inferred by demanding conservation of momentum. Let us define the fictitious 
inertial frame, I, in which M is stationary, and identify it with the frame in which 
M is actually stationary at the beginning of the scattering event (i.e. when m is 
at infinity). In this frame 

2M / 1 



X 



1 + 



(2.1) 



(Misner, Thorne and Wheeler 1973) where b is the impact parameter, and v re \ is 
the relative velocity in I. This relative velocity is given by 



w 



In manifestly covariant form 



— v| — |v x w| 
(1 - v.w) 2 

(V) 2 -1 



7 rel 



(w„v 



(2.2) 



(2.3) 



As w — > l,we find that v re \ — > 1; and in the limit that both w and v are much less 
than the speed of light, v re \ — > |w — v|. 

Now consider the scattering as viewed in the lab frame. We can write it as the 
product of a Lorentz boost, A, into I, followed by a rotation in 3-space by an angle 
X, and finally the inverse boost back to the lab frame, A -1 . The 4-momentum, , 
of the light particle is changed by an amount Aq^ 1 given symbolically by 



^^(A-^A-l)?' 



(2.4) 
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We shall use a convenient notation in which only two space dimensions are 
retained — one parallel to v, and one perpendicular to it. (By symmetry, only the 
first of these contributes to the dynamical friction, but the other is retained for 
clarity.) Thus we write 

( 1 \ 

= JwiTT- wcosO , (2-5) 
y w sin 6 J 

where 6 is the angle between v and w. For additional brevity, we define C = cos # 
and S = sin#. The boost A is thus written 

(7„(1 - vwC) 
j v (wC - v) | , (2.6) 
wS 

followed by the rotation 

(7„(1 - vwC) 
j v (wC — v)c — wSs I , (2-7) 
wSc + j v (wC — v)s 

where c = cos \ and s = sin%. The final boost back to the lab frame takes q^ 1 to 

(7^ (1 — vwC + v(wC — v)c) — j v vwSs 
-(I ((wC - v)c + v(l - wvCj) - j v wSs I . (2.8) 
wSc + j v (wC — v)s 

Thus we have, for \ small 

/ l 2 v v{wC - v)x 2 /2 - j v vwSx 
Aq' l =j w m\ jl(wC -v)x 2 /2- j v wS X \ (2.9) 
\ wS X 2 /2 + Jv(wC -v)x 

We now proceed to find the average rate of change of the momentum in the 
light particles, by averaging (2.9) over all possible impacts. If M were stationary, 
the number of collisions in the lab frame in a time interval dt, in the impact range 
[6, b + db] , would given by 

dN = nw2nbdbdt, (2.10) 

where n is the number density of light particles. To write dN in a Lorentz invariant 
way, first we note that in terms of the proper time interval, dto, in I, 

dt = j v dt . (2.11) 

Next, we examine the case that v is parallel to w, in which 

dN = j v j w (w — v)rio 2Trbdbdto } (2.12) 

where no is the proper number density of light particles in their rest frame. In 
terms of proper quantities and four vectors, and so in Lorentz invariant form, 
equation (2.12) may be written 
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dN = [{w^ f - l] 1/2 n 2Trbdbdt . (2.13) 
This may in turn be written in terms of lab frame quantities as 

dN = nV27rbdbdt, (2.14) 

where 

V ■ 



i I 2 I |2 

W — V — W X V 



1/2 



(2.15) 



Suppose that n is given by a distribution function /(w), such that the density of 
light particles is 

r -Y w mf(w)d 3 w, (2.16) 



then we find that 

^ = J 2TrbdbAq 11 Vf(w)d 3 w. (2.17) 

Integrating over b cancels the contributions proportional to \ by symmetry (the 
tedious inclusion of the third space dimension would make this explicit). Thus 



dq^ 



df J 47rAM 2 lw m I 7 l( w C-v) I y- + ^2~J Vf(w)d 3 w } (2.18) 

where A = ln(6 max /6 m i n ). Since we are only considering weak scattering, A is 
independent of w. For a discussion of a natural choice for 6 m i n in the case that 
M is an extended distribution of mass see White (1976). We now concentrate 
on equation (2.18) in two interesting limits: first we re-derive the Chandrasekhar 
formula; and then we find the ultra-relativistic limit as w — > 1 for general v. 

In the limit that w and v are much less than the speed of light, the 3-space 
part of (2.18) is equivalent to 

^ = / ^AM 2 m W ~ V 3 /(w)t/ 3 w. 
at J |w — v| 

For an isotropic distribution /, 



dq ^ AM 1 raw 

dt v 3 



/(w)d 3 w, (2.19) 



w<.v 



which is the usual Chandrasekhar (1943) result. 

In the limit that w — > 1, we must also demand that m — > 0. We write / in 
terms of q, by letting j w m — > q } whence 

da* r (il v (c-v) 

q 1 IQttAM 2 ll(C-v) I (1 - U C)g/(q)(/ 3 q. 



dt j x S 



For an isotropic distribution 
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% = J^M* P1 tv \l |. (2.20) 

The momentum, of M will suffer an equal and opposite change, so 

dQ» _ dq» 
dt dt 

From this we may immediately find the interesting limit of the friction formula in 
which M is a hard photon. Taking the limit v 
we see that 

= T* AQ p 

(which is very small in any astrophysical case). 
Noting that 

Q» = lv M | v | , (2.23) 




we derive 




it < u n 1 1 1 • < 2 ' 24 > 

and thus that the energy and momentum parts of (2.18) are consistent, given that 
the component perpendicular to v must vanish by symmetry. 

Thus finally, in the case w — > 1 and where / is isotropic we have, from (2.20) 
and (2.24), 

dv 64 Mo , 

— = ttA— v. (2.25) 

at 3 7„ 

Re-inserting the dimensional constants G and c this is 



dv 64 t G 2 Mp 

— = TTA 

dt 3 7„c 



ttA _ 3 r v. (2.26) 



Comparing equations (2.25) and (1.2) we see that, (i), the fact that all the light 
particles have higher velocities than M does not reduce the dynamical friction; 
and (ii), the rate of change of v due to dynamical friction is modified by a factor 
of roughly 16/37„ over the simple non-relativistic case. 



D. Syer: Relativistic Dynamical Friction 



6 



3 Discussion 



We should be aware of the effect of direct collisions between M and the light par- 
ticles. If the light particles are photons, this corresponds to the radiation pressure 
force on M, assuming it is optically thick. Suppose that M was a 'shiny disc' of 
radius a, meaning that, in a direct collision, the light particles are elastically re- 
flected in the direction of v. Going through a similar derivation to the above we 
would find that, in the limit w — > 1, 

dv 16 p 

~r = Ti" — 77V. 3.1 

dt 3 7„M v ; 

The ratio of the dynamical friction rate to the direct friction rate is A AM 2 j a 2 ~ 
ylUg SC , where v esc = (2M / a) 1 ! 2 is the escape speed from the surface of M. 

For objects, such as stars, v esc <C 1 and hence dynamical friction by photons 
can never be important. There are, however, at least two cases in which dynamical 
friction against relativistic light particles could be interesting: if v esc ~ 1; or if the 
optical depth, r, of M to the light particles satisfied r <C ^esc- 

If v esc ~ 1 then the weak-scattering approximation would no longer be valid, 
but actually the dynamical friction would be greater (the effective A would be 
vary large). Thus neutron stars and black holes would feel a greater force due to 
dynamical friction than due to radiation pressure, if immersed in a background of 
photons of lower energy. The ratio of the radiation pressure to the gravitational 
force on a black hole of mass M in the vicinity of an object radiating at its 
Eddington limit is 

-^radiation r 1n -l5 M fo o\ 

— L M ~ 10 — - , (3.2) 

r gravity 

where Lm is the Eddington luminosity of M itself in natural units. Thus for 
radiation pressure to have an appreciable effect on the orbit of M in a Hubble 
time, we require that 

1Q _ 15 M tHubble ^ (3 3) 

Mq torbit 

where t rbit is the period of the orbit. Writing t rbit in terms of the separation r/M 
of two equal mass black holes, the inequality (3.3) becomes 

V 3/2 £ 10 7 , (3.4) 

independent of M. For r/ ~ 1 gravitational radiation would take over as the domi- 
nant effect on the orbit, but gravitational radiation decreases quickly with increas- 
ing separation, so radiation pressure might be important on timescales as short as 
10 6 years. Dynamical friction would be comparable or larger. 

If the light particles were neutrinos, then a star would be very optically thin 
and direct collisions would not provide any friction at all. Neutrino balls (Holdom 
1987, Holdom and Malaney 1994) are an example in which the neutrino density 
would be high enough to affect the orbit of a star. In such objects the neutrinos 
are not strictly massless, but they are degenerate, and so they are virtually all 
ultra-relativistic. 
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